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.^ , We pose a novel cellular network planning problem, considering the use of renewable energy sources 

^^ . and a fundamentally new concept of energy balancing, and propose a novel algorithm to solve it. In terms 

of the network capital and operational expenditure, we conclude that savings can be made by enriching 

cellular infrastructure with energy harvesting sources, in comparison to traditional deployment methods. 
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,-j-' ' Although equipping base stations (BSs) with Renewable Energy Sources (RESs), such as solar panels 

^r) • and wind turbines (either to support or replace traditional electric grid connection with energy harvest- 
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ing [|2]]), has been proved technically feasible, the research on RES-powered cellular networks is still 
in its infancy. One of the fundamental questions to be answered is how to plan the topology of RES- 
enabled cellular networks (i.e. physical location of BS and the location of RES) to satisfy coverage 
in accordance with quality of service (QoS) needs (defined in some specific sense), while significantly 
reducing capital and operational expenditure (CAPEX and OPEX, respectively). Therefore, inspired by 
the studies of BSs deployment in traditional cellular networks Q, |@], this letter attempts to answer this 
question by formulating a novel optimization framework for planning the RES -enabled cellular network. 
As the problem is shown to be NP-hard, we build a heuristic for cellular topology planning consisting 
of two phases: (i) QoS-aware BS deployment and (ii) energy balancing connection. Numerical results 
show that the proposed heuristic brings CAPEX and OPEX savings in comparison to traditional network 
deployment methods. 

II. System Model 

Network Characterization: We consider an area where a cellular network must be deployed or 
enhanced/upgraded with new features. As downlink is the bottleneck of a cellular network (since far 
more traffic is sent over downlink than uplink) it needs to be handled first in network planning. We thus 
focus on downlink transmission. Then we assume a set of candidate sites for BS installment labeled 
M = {1, ...,N}, with a respective installation cost c n , n G J\f. Further, we assume a set of test 
points (TPs), i.e., centroids, where a given amount of traffic or a certain level of QoS service must be 
guaranteed (3l Sec. Ill- A], located at sites M. = {1, . . . ,M}. Expressly, TP can be interpreted as a (set 
of) user equipment in the cellular networks (with certain properties, e.g. minimum required signal to 
noise ratio, defined individually per each TP). 

Energy Consideration: As BS n is connected to an electric grid consuming power Pn , Mn £ J\f, 
we aim to minimize the energy consumption at the network planning stage. We thus assume that each BS 
is equipped with RES, e.g., solar panel, wind turbine. The instantaneous (and aggregated from all energy 
sources) power delivery capacity of RES for BS n, Z n , is a stationary stochastic process described by 
PDF fn (z n ). As a fundamental novelty, for potential energy delivery improvements, we assume that 
BSs can have dedicated power line connections for energy balancing^ from other RESs. On the other 
hand RESs are not equipped with any battery to store energy (as fixed electric grid connection guarantees 
a constant energy supply in case of insufficient energy from RES). This assumption results in a lower 
bound for the system planning, noting however that our model can be extended to consider energy storage. 

'Note: this concept has been independently presented in (5) once our paper was already under review. 



The power transferred from BS t to BS n (through its co-located RES) is denoted as Pt tU , while the 
associated (distance-dependent) cost, e.g. installation, maintenance, operation, of connecting BS t and n 
is denoted as ct >n , t G J\f. 

Network Planning Constraints: We define the following Boolean variables. The BS deployment 
indicator b x = 1 if a BS is installed on site x G A/", and b x = otherwise. The TP assignment indicator 
Pm,n = 1 if TP m € .M is served by BS n G A/", and p m;n = otherwise. The inter-RES power line 
connection indicator for each BS Ct n = 1 if BS £ G .A/", n G TV are connected by a power line, and 
Q,n = otherwise. We then have 

Pm,n < b n ,Cj jn < b n ,Cj jn < b t . (1) 

We further assume the following constraints. Each TP has to be connected to exactly one BS, i.e., 

^p mi „ = l,VmeM. (2) 

neJV 

The number of TPs served by each BS is upper bounded, such that 

Y, P m ,n<£b n ,VnG/V, (3) 

where B is, without a loss of generality, the maximum number of TPs served by each BS, i.e. total radio 
resources available to each BS. Each BS is assumed to occupy the entire spectrum, i.e., factor- 1 reuse 
strategy is adopted. 

Propagation Environment: We use H = \hm,n{Mx.N to denote the channel-dependent propagation 
gain matrix where h m>n is the long-term propagation gain of the link between TP m (averaged over all 
users in TP m) and BS n [3j Sec. Ill- A]. Denoting 5 m as the thermal noise power at TP m and P n as 
the total transmission power of BS n for each associated TP, we define SINR between TP m and BS n 
as SINR m , n = P n h m , n \T, t ^ n ptPthm,t + 5 m ) , V?n G M, n G M, where p t = ^ jeM Pj,t/B denotes 
the ratio of occupied resources to all available resources at BS t (5J Eq. (10)]. In addition, TP m has a 
minimum QoS requirement defined in SINR terms, 7 m , such that SINR mn > p m ,n7m> Vm G M, n G N. 
For the ease of subsequent discussion we transform the non-linear QoS constraint into an equivalent 
linear formulation through Sec. IV-C] as 

SINR m , n > p m ,nlm «=> 

-'"r?i,n(-t Pm,nJ ~r J~^n'T"m,nPm,n ^ 

| B^ 1 J2 Yl Pj,t P thm,t + Sm , (4) 

t^njeM J 



where M mn = 7 m ( Yltj^n Pth m ,t + <^m ) is a sufficiently large constant such that d4) holds for any 
Pm,n G {0, 1}, V?n G M, IJ£ .A/". 

Power Outage Constraints: The energy consumption of each BS is limited by the sum of its all 
energy sources. Thus, an energy outage happens at BS n when its energy supply is insufficient for the 
energy demand to serve its users. Setting a very small upper bound of the energy outage probability for 
each BS n, (p n , (to limit energy supply loss due to, e.g., blackout) yields the following constraint 

Un I / j Pm,n-fri. t ±. n r, n -\- > J~^n,t^-n,t 

J](l - Et,n)Pt,nCt,n ~ ff (<Pn) ) < 0, Vn G M, (5) 



where e in G [0, 1] is the energy loss factor on the power line between RES t and no, -P« is the 
required static power for operating BS, e.g., cooling, baseband processing (H, and In {■) is the inverse 

of 4 p) oo = /*•/£> m. 

III. Problem Formulation via Optimization 

The objective is to select a subset of candidate BSs within J\f and to assign TPs within M. to an 
available BSs taking into account the QoS requirement, BSs installation, inter-RES connection cost, 
and the electric grid consumption. Combining the constraints from Section JI] we obtain the following 
optimization problem 

min } c n ,b„ + } ct, n ct,n + 

{p m ,„,c t ,„,b„,P^ ) >0,P t ,„>0} f~ 




XT J] P n )b n + 2 £ t,nPt,nCt,n (6a) 

\n£j\f t,n£j\f J 

V v ' 

(iii) 

subject to O, ©, ©, ©, ©. (6b) 

The terms (i) and (ii) correspond to the total installation and the connection cost, respectively. The 
term (iii) characterizes the cost of consumed power from the electric grid, where A is the energy price 
expressed in cost unit/kWh, and T is the specified life cycle of the deployed cellular network. Sum of 

2 With Et, n we guarantee that we do not deploy inter-RES links to transfer negligible amounts of energy between RESs. 



(i)-(iii) terms in (I6ab denotes CAPEX and OPEX. Note that consideration of the random energy arrivals 
to the harvesters in our planning framework are considered in (l6bl through ©. 
Proposition 1: Problem © is NP-hard. 

Proof: (Sketch) Following Q, the problem of © can be reduced to an uncapacitated facility location 
problem which is NP-hard [3, Eq. (3)-(6) and Sec. IV] and implies NP-hardness of ©. □ 

IV. Problem Solution via Heuristic 

To solve ©, we propose a heuristic that decomposes © into two decoupled subproblems. The first 
subproblem is the QoS-aware BS deployment (that deploys BSs and connect TPs to BSs according to the 
energy distribution of BSs, assuming all candidate BSs are disconnected). The second subproblem is the 
energy balancing (that further reduces the total cost of network operators by balancing the benefits given 
through energy sharing among deployed BSs and their RESs and the incurred cost on BSs connection). 

A. QoS-Aware BS Deployment 

First, we assume that no inter-RES connections are possible, i.e., c t ,n = Pt,n = 0, Vn, t G Af. Defining 

Ci(v, w) = ]C neA f c n v n + XT J2 n( z_\r w n y n , where v = [v n ]jvxi, w = [w n ] Nxl , © reduces to 

min Ci(b,P) (7 a) 

{p„ l ,„,b„,P,<' ) >0} 

subject to CD, ©, ©, ©, (7b) 

b„ ( Yl Prn,nPn + P<?> ~ P^ ~ ff (<Pn))< 0, (7c) 

\m£M ) 

where b = [b n ]jv x i, P - [Pn^Nxi- 

Proposition 2: Suppose an optimal solution to ([7]) exists, denoted as a 3-tuple {b*,p*,P*}, where 
P = [p m ,n]MxN and * symbol denotes optimality. Defining A n (u) = - Y /meM u m>n P n -Pi' ] +I^\ip n ), 
where u = [u min ] MxN , then 

P^)*b*=ma X {-A n (p*)b;,0} ; (8) 

where A n (p*) is the remaining energy from RES applied to transmission and operation of BS. 

Proof: For b* = 1, (TTcl ) reduces to P„ > — A n (p), which together with P„ > implies that 

P^ > max {-A n (p), 0} , Vn € M. (9) 



As the objective function of ([7]) is strictly increasing in PJf , Vn G M, when b* = 1, we conclude 
that ([9]) will be always tight in the optimal condition, i.e., Pn = max{— A n (p*), 0}, \fn € TV. On the 
other hand, © holds when b* = 0. D 

Applying © to © yields 

nun } XT £ max{|^ - A n (p), |^} b n (10a) 

subject to d7b]l. (10b) 

Notice that solution to (fTOl is equal to the solution to the original problem ©. From (TTOl follows that it 
is beneficial to deploy BSs at sites with positive A n (p), and to connect TPs to the deployed BSs with 
surplus renewable energy. We devise a heuristic to solve (ITOl (i.e. an algorithm that reaches as close 
as possible to a final solution) which is presented in Algorithm [Q Therein, (for ease of exposition) the 
total radio resources are assumed sufficient for network planning, i.e. BN > M. Note that Algorithm 1 
might not output the assignment for all existing TPs when, e.g., QoS requirements are too strict. In this 
case, e.g., new positions (or the number) of BS must be considered and Algorithm 1 must run again. The 
convergence for Algorithm [Q can be achieved in at most N outer loops, where the worst case represents 
reaching M = 0. 

B. Energy Balancing Inter-RES Connection 

Let M denote the set of deployed BSs based on b, p, P (i.e. the output of Algorithm [T). Defining 

C 2 (x,y,z) = XT ^ tn£ tf (et,nZt,nCt,n + ^C tt n)+Y,n&Af (Vn+w)}' Wnere X = [ x t,n]NxN, Y = 

[yn]Nxi, and z = [zt t n]NxN- After the QoS-aware BS deployment shown in Section HV-Al © further 
reduces to 

min C 2 (C,P,R) (11a) 

{c t ,„,P^ 9) >0,P t ,„>0} 

subject to © with p„ v „ from p, replacing M by J\f, (lib) 

where C = [c ttn ] NxN and R = [Pt^NxN- 

We propose a second heuristic (algorithmic solution) given in Algorithm |2] to reach as close as possible 
to a solution of (TTTT t. The heuristic starts with a fully connected RESs topology, with a set of inter-RES 
connections denoted as £ = { e t,n} tn( zj^- Within each iteration of the heuristic, the inter-RES power 
lines transferring the least amount of energy are removed from £ until no cost saving can be achieved. 
Notice that {C = [l]vt jn ;P;R = [0]vt,n} is a feasible solution to (fTTT ). thus the feasibility of Step [2] in 



Algorithm 1: QoS-aware BS deployment 



Input: \/m,n,t: M; M\ ct,„, P n , S m , ip n , c„, Pi°\ B (s.t. BN > M), e t , n , lm\ fi r) (z n ); H 

1: connect all TPs to closest BSs with free (unassigned) radio resources 

2: fi := {m\ X^ neJ v-SINR m ,„p m ,„ < j m ,m G M}, i.e. infeasible TPs 

3: if n = then 

4: F:=l;p:=p;b:=b;P:=P 

5: repeat 

6: from BS set with least number of TPs, i.e., n G argmin nG Af X^mewi Pra,m randomly select one and disconnect its 

associated TPs, i.e., M 

7: for all TP m G M<, do 

8: find b n by assigning TP m to the closest BS n (except n ) with positive A n (p) without violating $3$ 

9: if such b n exists then 

10: p m , n = 1, n 6 AT; M, := 7W \{TP m} 

1 1 : else 

12: p m ,„ = 0, n G Af 

13: end if 

14: end for 

15: n := {m| I] nGAr SINR m>n p m , n < 7 m , m G A4} 

16: if n = and X = then 

17: Af := JJ\ {no}; p := p; b := b; P := P 

18: end if 

19: until Q / or J\f = 

20: else 
21: F:=0 

22: end if 
Output: F := 0: no solution found; F :— 1: p := p, b := b, P := P, i.e. optimized BS configuration 

Algorithm |2] always holds. In addition, termination criteria for the loop (Step [T]-[5]> in Algorithm |2] are 
achievable within \£\ loops. 

C. Heuristics Complexity Analysis 

The complexity of Algorithm [T] is determined by the double-nested loop (Step|5rfT9l) whose complexity 
is 0(MN 2 ) = 0(N 3 ) since M = 0(N), i.e., the number of TPs is limited by M < BN. The complexity 
of Algorithm |2] is governed by the time needed to solve a Linear Program (LP) (StepE). As each LP can 
be solved in polynomial time [9, Theorem 8.5] its complexity is 0(X t ), where X is the number of LP 
variables and t characterizes the running time of the LP solver. Thus the complexity of Algorithm 12 is 
0{X)0{X t ) for X = N 2 . In turn, the complexity of the complete heuristic (Algorithm [J and Algorithm^ 
is 0{N 2t+2 ). 



Algorithm 2: Energy Balancing inter-RES Connection Deployment 



Input: complete set of inter-RES connections £ ; Vt,n: et, n , c„; p,b,P from Algorithm [2 § = oo 

1: repeat 

2: find C, P, R, i.e. the solution of Jilt , with p, b, ct, n = 1 Ve t , n £ £ 

3: target zero and the smallest positive Pt. n in R, then remove corresponding connections from £ , i.e., £ :— £\{et,n} 

4: §:= min{§, C 2 (C,P,R)j 

5: until C 2 (C, P, R) > § or £ = 

6: if Ci(b,P) <Sthen 

7: £ := 

8: end if 
Output: optimized inter-RES connections according to £ 

V. Numerical Results 

As an illustrative simple example we consider a 3 km x 3 km area where m TPs are arbitrarily located 
and n BS candidate sites (forming a grid) are identified to provide services for TPs. For a cost unit 
of € we assume all RES-enabled BSs have fixed and equal installation cost, c n iflOl Sec. II-B] and 
fixed transmission power ||7] Sec. II-C]. Furthermore, for a fair comparison we assume that all traditional 
(no-RES enabled) BSs have an installation cost of c n and find optimal cellular network structure using 
Algorithm [T] considering no RES and no inter-RES connections. Assuming distance-dependent pathloss- 

-{L A +L B losi (d m ,„/km)} 

only scenario, the elements of H are given as h m>n = 10 "> [61 Eqs. (1), (7)], where 

d m>n denotes the distance between TP m and BS n and La and Lb are the empirical constants provided 
by El. 

As the distribution of harvested power at location n is scenario-dependent, for simplicity we assume 
it being uniformly distributed, i.e., fr (z n ) = , jj; , \/z n € [a n ,b n ], n G A/", with a n and b n being the 
minimum and maximum harvested power at location n, respectively (note that we can use any continuous 
energy distribution in our model). All deployment parameters are summarized in Table [Q representing 
values describing typical network scenario, following, e.g. O, H, ifTTl . Numerical results are generated 
based on the method of batch means with 100 simulation runs for the confidence level of 95%. Due to 
space constraints we focus on two most representative cases. 

A. CAPEX and OPEX versus Network Life Cycle 

The result is presented in Fig. |l(a)| For each simulation run new location points for BS and TPs (with the 
respective H) have been randomly generated. We immediately observe that the introduction of RESs into 
network planning brings cost saving to the operator, compared to the traditional deployment structure. 
This benefit is slightly boosted by the energy balancing among RESs, as for the considered network 



TABLE I 

Parameter Values Used in the Numerical Evaluation 



No. candidate sites 


N = 9 


Network life cycle 


T = [6, 20] years 


No. TPs 


M = 20 


Electricity price 


A = 0.3€/kWh 


No RES-based BS depl. 


c n = 55 k€ 


RES-based BS depl. 


c„ = 60 k€ 


Inter-RES conns, cost 


<H,n = 10€/m 


BSs oper. power 


P^ o) = 19 W 


Transmission power 


P„ = 20W 


No. res. blocks 


B = {6, 12} 


Thermal noise power 


S m = -114dBm 


SINR requirement 


7m = 0dB 


Propagation coeff. A 


L A = 148.1 dB 


Min harv. power 


a n = [0, 100] W 


Propagation coeff. B 


L B = 37.6dB 


Max harv. power 


6„ = [100, 200] W 


Power outage prob. 


fn = 5% 


Energy loss 


£t,n = 1% 



coeff. -coefficient; conns. -connections; depl. -deployment; harv.-harvested; no. -number; oper. -operational; prob.-probability; res.-resource; 

reported values are equal Vm, n, t 
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(a) x: No REs (optimal depl. via FICO®Xpress); o: No inter- (b) x: No inter-RES conns. (B = 6); o: Inter-RES conns (B = 6); 
RES conns. (Alg. [TJ; *: No RES (Alg. [[}; o: Inter-RES conns. *: No inter-RES conns. (B = 12); o: Inter-RES conns. (B = 12) 
(Alg[j]+Alg.r2) 

Fig. 1. CAPEX and OPEX as a function of: (a) network life cycle T, for randomly chosen a n , b n from Table U for each 
location n and B = 12; (b) harvested energy spread max6„ — mina n for T = 10 years, where for each simulation point (from 
left to right) k — 0, ... ,9 : E[z n ] = 145 — 5fc (W); whiskers denote 95% confidence interval; alg.-algorithm; depl.-deployment; 
conns.-connections. 

configuration, per BS, most of the energy is provided by the RES alone. Notice that CAPEX and OPEX 

with inter-RESs connections does not vary with the network life cycle, as the whole network is basically 

self-powered (without grid power). Comparing the effectiveness of the developed algorithmic solution 

(considering the lack of RESs) with the solution to ©, computed using FICO® Xpress Optimization 

Suite version 1.23.00, for 10 independent realizations of H, we conclude that our heuristic is close to 

the optimal deployment. 
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B. CAPEX and OPEX versus Harvested Energy Spread 

The result is given in Fig. |l(b)| for one randomly generated set of TP locations, with the respectively 
generated H (one time only). CAPEX and OPEX roughly increases as the difference between mina n 
and max b n becomes larger, irrespective of B. This is because of strong randomness in available energy, 
which in turn incurs extra cost on electricity supply (no inter-RES connections) or power lines (inter-RES 
connections). In addition, we also observe that deploying inter-RES connections with large B is the most 
cost-effective option for network operators. This is because BSs with large B normally allow the small 
number of deployed BSs thus saving the deployment cost, while also lead to large energy deficiency at 
some 'crowed' BSs. This is especially visible with high variance of harvested power, which increases a 
benefit of RES use by deploying inter-RES connections. 

VI. Conclusions 

In this letter we have developed a novel cellular network planning framework considering the use 
of renewable energy sources and energy balancing. For the posed problem we have developed a novel 
heuristic. Our numerical results demonstrate CAPEX and OPEX savings in comparison to traditional 
deployment strategies. 
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